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ABSTRACT 



Aims. We study the content in S-bearing molecules of protoplanetary disks around low-mass stars. 

Methods. We used the new IRAM 30-m receiver EMIR to perform simultaneous observations of the 1 lo - loi line of H2S at 168.8 GHz 
and 223 - I12 line of SO at 99.3 GHz. We compared the observational results with predictions coming from the astrochemical code 
NAUTILUS, which has been adapted to protoplanetary disks. The data were analyzed together with existing CS J=3-2 observations. 
Results. We fail to detect the SO and H2S lines, although CS is detected in LkCalS, DMTau, and GOTau but not in MWC480. 
However, our new upper limits are significantly better than previous ones and allow us to put some interesting constraints on the 
sulfur chemistry. 

Conclusions. Our best modeling of disks is obtained for a C/0 ratio of 1.2, starting from initial cloud conditions of H density of 
2 X 10^ cm"' and age of 10^ yr The results agree with the CS data and are compatible with the SO upper limits, but fail to reproduce 
the H2S upper limits. The predicted H2S column densities are too high by at least one order of magnitude. H2S may remain locked 
onto grain surfaces and react with other species, thereby preventing the desorption of H2S. 



Key words. Stars: circumstellar matter 
Submillimeter:planetary systems 



planetary systems: protoplanetary disks - individual: MWC480, LkCal5, DM Tau 



1. Introduction 

Understanding the evolution of gas and dust accretion disks 
around young stars is one of the prerequisites for handling the 
processes leading to planet formati on. The temperature struc- 
ture of disks is roughly und erstood (Chiang & Goldreichl l 19971: 
iMen'shchikov & Henningll997 : D'Alessio et al. 1999.) . but very 
fundamental disk properties, such as the ionization degree and 
the dust and gas mass distribution, are poorly constrained. As 
no direct tracer exists, the derivation of the gas mass distribution 
requires in-depth knowledge of the disk structure and chemistry. 
Therefore, studying the chemistry becomes a necessary step, be- 
cause molecular abundances and gas distribution are strongly 
coupled in both the modeling and analysis of observational data. 

Vertically, the current paradigm of the chemistry of 
protoplanetary disks is a layered model (see for example 
Ivan Zadelhoff et al1l2001h . In the upper layer, the large incident 
UV flux from the central object results in a photodissociation- 
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dominated layer (PDR) with a chemistry in equilibrium. The 
PDR size evolves wit h grain growth, since UV flux can penetrate 
deeper into the disk ( Aikawa & Nomurall200"6l : IChapillon et al.l 
I2OO8I: IVasvuninetal.l l201l'). Below the PDR, the molecules 
are expected to concentrate in a warm molecular layer. All 
chemi cal models published so far have produced si milar re- 
sults dAikawa & Nomural I2OO6I: ISemenov et alJl2005h . Finally, 
because of the large extinction, the outer disk midplane is cold 
(~ 10 K) and molecules are expected to stick onto dust grains. 
The chemistry is then rather similar to what is encountered in 
cold dense cores. Radially the "snowline", which is defined by 
the radius where H2O molecules start to evaporate from the ice 
mantle of dust grains, delineates the zone where rocky plan- 
ets can be formed. For a TTauri star similar to the young Sun, 
the "snowline" at midplane should typically be located around 
-0.5-2 AU. 

The ideal scenario presented above is loosely constrained. 
There are already several observational facts that reveal deficits 
in current chemical models for the outer disk that contains 
the reservoir of gas and dust mass (~ R> 30 AU). Very re- 
cently, sensitive observations with Herschel/HIFI have revealed 
that the emission o f H2O is significantly weaker than predicted 
(iBergin et alJl201(il) . Another example is given by our mm ob- 
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servations, which show " cold" molecules (at temperature lowe r 
than ~ 10 K) such as CO (iDartois et alJl2003l: lPietu et al. ll2007h 
C2H (|Henningetai1 l20"l0h . or CN and HCN (Chapill on et al.l 
1201 Ih . In the last cases, the disk density exceeds the critical den- 
sity of the o bserved transitions and subth ermal excitation cannot 
be invoked jPavlvuchenkov et al.ll2007l) . The usual explanation 
given is the role of turbulence. In theory, the vertical and radial 
mixing in the disk s hould allow a partial replenishment of the 
cold midplane layer (ISemenov et al. I I2OO6I). However, this does 
not appear to be sufficient by itself ( Hersant et alJ l200"9l) . The 
turbulent mix ing efficiency actual ly depends on the photodes- 
orption rates dHersant e t al."2009') which may be much higher 
than originally thought (iOberg et al. .2009) . 

For inner disks, where planets should form (~ R< 30 AU), 
some very interesting unresolved observations (mostly coming 
from the Spitzer satellite) begin to unveil the chemical content of 
planet-forming regions. These observations show that not only 
CO but also many other molecules such as H2O, HCN, CO2 , 
C2H2, OH, etc. . . are present in inner disks dCarr & Naiital20()8h . 
So far, there has been no H2S detection. Interestingly, a recent 
paper from[Glassgold et al. ( 2009) suggests that there is no need 
for transport in the inner nebula to explain the presence of water 
vapor at the surface of the very inner disk. The disk irradiation by 
UV and X-rays may be sufficient to create water "in situ", in the 
warm disk atmosphere. This model accounts rou ghly for the IR 
water vapo r emission observe d in TTauri disks bv lCarr & Naiital 
(l2008h and lSalvketai] (l2008h without radial or vertical mixing 
in disks. 

R ecent models of so lar nebula investigating the sulfur chem- 
istry dPasek et all l2005h within the 20 central AU found that 
the main gas component of sulfur-bearing molecules would be 
H2S, while it should be in the form of FeS, MgS, or CaS in 
solid state (meteorites). The abundance of those species are ex- 
pected to be variable with respect to the time evolution of the 
snowline and the water vapor content. Moreover, sulfur-bearing 
molecules such as H2S and CS (and SO to a lesser extent) are 
relatively abundant in comets. For example, in Hale-Bopp, the 
abundance of HtS is 10% that of CO. H2S is supposed to come 
from the evaporation of ices, wh ile CS and SO are the daughter s 
of CS2 and SO2, respectively dBockelee-Morvan et a n 12004 . 
Nonetheless, the only S-bearing molecule that is firmly detected 
in cold outer TTauri disks is CS. 

Several transitions of CS, the most abundant S-bearing 
species i n cold molecular cor es (Dickens et al. 2000), were de- 
tected bv iDutrev et alj dl997h . who found an abundance (~ 3 ■ 
10 '°), about 30 times lower than in TMCl. Such a difference 
is not surprising since the gas and dust have evolved inside the 
disk compared to dense cores. 

In this paper, we go one step further by studying other S- 
bearing molecules using the IRAM 30-m radiotelescope. The 
observations and results are described in S.2, the chemical mod- 
eling is presented in S.3, and we discuss in S.4 the implications 
of our non detections before we conclude in S.5. 



2. Observations and results 

2.1. Sample and observations 

Many molecular lines have been observed around MWC 480, 
LkCalS, and DM Tau dDutrev et al.lll997h . and the disk struc- 
tures are well known dPietu et al.ll2007h . DM Tau and LkCalS 
are TTauri stars of 0.5 and 1.0 Mq, respectively, while MWC480 
is an He rbig Ae star of 1.8 M p. The CO disk of GO Tau is re- 
ported in lSchaefer et alJd2009l) . It appears to be a good candidate 



Table 2. Physical parameters used to derive the best fit models 
and the surface density upper limits 



Source 


MWC 480 


LkCa 15 


DM Tau 


GO Tau 


inclination (") 


38 


52 


-32 


51 


P.A.C) 


57 


150 


65 


112 


V,„,(km.s-') 


5.10 


6.30 


6.04 


4.89 


Vioo(km.s"') 


4.03 


3.00 


2.16 


2.05 


5,,(km.s-') 


0.2 


0.2 


0.2 


0.2 


Tioo(K) 


30 


15 


15 


15 


R„„(AU) 


1 


45 


1 


1 


R«„(AU) 


500 


550 


800 


900 


CS R<„„(AU) 






560± 10 




CS p 






0.75+0.05 


1.0±0.2 



Notes. P.A. is the position angle of the disk rotation axis, ; the incli- 
nation angle of the disk (0 means face -on), Wsvst the systemic velocity. 
The velocity laws (Vioo) are Keplerian jPietu et al.l2007l : [Schaefer et all 
i200g). 6„ is the turbulent line-width component. The temperatures 
and surface densities follow the laws T(r) = rioo('7100AU)"''^ and 
S(j-) = E3oo(f/300AU)"' \ respectively. For CS, in the case of DM Tau, 
the outer radius and the surface density radial dependence p derived 
from the data have been used. This is also the case for p in GO Tau. 
Derived parameters are given with their errorbars (a "-" means that the 
standard values are assumed.) 



for a deep search of chemistry since it exhibits a large disk (CO 
outer radius R,,,,, ^ 900 AU) that is very similar to that of DM 
Tau. Table[T]summarizes the stellar properties of the four stars. 

The observations were performed using the new 30-m het- 
erodyne receivers (EMIR) mid-May 2010 under good weather 
conditions. We used the wobbler switch mode with offsets of 
60" since the sources are unresolved in the telescope beam, pro- 
viding spatially integrated spectra. We observed the H2S 1 lo-loi 
line at 168.763 GHz simultaneously with the SO 223 - I12 line 
at 99.299 GHz using the correlator VESPA with a spectral reso- 
lution of 20 KHz. Venus was regularly checked for pointing and 
focus corrections. The forward efficiency was set to 0.95 and 
0.93 at 99 and 168 GHz, respectively. We used flux density con- 
version factors (Jy/T*) of 6 Jy/K at 99 GHz and 6.5 Jy/K at 168 
GHz. We took the values recommended by the IRAM staff (see 
http ://w w w. iram .es/IRAMES/ main Wiki/Iram3 OmEfficienc ie s ) . 

The CS J=3-2 data of GO Tau were observed in 2009 at 
the IRAM 30-m radiotelescope using EMIR and reduced as de- 
scribed above. The CS J=3-2 data of MWC 480 and LkCal5 
were obtained with the IRAM 30-m telescope in 2000 si- 
multaneously w i th the DCO^ J=2-l observations published in 
iGuilloteau et alj d2006h. usir ig the same observational strategy 
as described in lDutrev et"aD dl997h . For DM Tau, CS J=3-2 was 
mapped with the IRAM array in the course of an imaging project 
(Guilloteau et al., 2011, in prep.) and we only present here the 
integrated spectrum (integrated over the CO outer radius). 

All the 30-m data were reduced using CLASS. The DM Tau 
CS map was analyzed using CLIC and MAPPING and an inte- 
grated spectrum was produced. 



2.2. Column density derivation and results 

Figure [T| shows the integrated spectra obtained for all the 
sources. Neither H2S nor SO are detected but the observations 
provide for the first time meaningful upper limits, significantly 
better (by a factor ~ 7 - 10) than those quoted in JJutrev et afl 
(I2000h . CS J=3-2 is detected in three disks: GO Tau, LkCal5, 
and DM Tau, but not in MWC 480. 
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Table 1. Sample of stars and stellar properties. 



Source 


Right Ascension 


Declination 


Spect.Type 


Effective Temp. 


Stellar Lum. Stellar Mass 


Age 


UV flux 




{",'"), (J2000.0) 


(*,'",') (J2000.0) 




(K) 


(Lo) (Mo) 


(Myr) 


(Xo) 


LkCa 15 


04:39:17.78 


22:21:03.34 


K5 


4350 


0.74 1.01 ±0.02 


3-5 


2550 


DMTau 


04:33:48.73 


18:10:09.89 


Ml 


3720 


0.25 0.53 ± 0.03 


5 


410 


MWC 480 


04:58:46.26 


29:50:36.87 


A4 


8460 


11.5 1.83 ±0.05 


7 


8500 


GOTau 


04:43:03.05 


25:20:18.8 


MO 


3850 


0.37 0.5 


3 


410 



Notes. Coordinates J 2000.0 deduced from the fit of the 1.3mm continuum map of the PdBI performed in'Pietu et alj ( l2006h . Colum ns 3,4, 5, 6, 
and 7 are the spectral type, effective temperature, stellar luminosity, and age as given in^Simon et al. (2000). Masses are ta ken from lPietu et al.l 
( I2007I) . for DM Tau, MWC480, and LkCal5. The stellar UV luminosities that are given in C ol. 8 in units of the Draini | |1978|) interstellar- UV field 
are taken from lBergin et all l l2004h (LkCa 15 and DM Tau) or computed from the lKuruczl i 19931) ATLAS9 of stellar- spectra (MWC 480 and GO 
Tau). They are given for a distance of 100 AU from the star. Half of this flux is supposed to scatter downwards inside the disk. 

LkCa 15 MWC 480 DM Tau GO Tau 




-10 10 

Velocity (km/s) 



Fig. 1. Observations of SO 223 - I12, H2S lio - loi and CS 3-2 in the four disks. For CS J=3-2, the best models have been 
superimposed for all sources. For SO and H2S, the models always correspond to the 3cr upper limits. 



Deriving proper estimates of upper limits is complicated in 
the case of Keplerian disks because of the existence of veloc- 
ity, temperature, a nd density gradients dGuilloteau et al.ll2006t 
iDutrev et alj|2007l) . For chemistry, since the surface density and 
the temperature have radial dependencies, defining an averaged 
column density over the whole disk implies to merge chemically 
different regions. This can only be properly done by taking the 
existing gradients into account. Because the characteristic size 
of the disks is several 100 AU, and is unresolved in the 30- 
m beam, our observations are most sensitive to the 200 - 400 
AU ra nge, depending on the radial distribution (see iPietu et al] 
l2007h . We thus chose to proceed as follows. We used DISKFIT 
(a dedicated radi ative transfer code to protoplanetary disks, see 
iPietu et ai]|2007h to generate integrated spectra and derived the 
best model by adjusting disk parameters. The models provided 



by DISKFIT were compared to th e observed spectra using the 
minimization routine described in iPietu et al] ( l2007h . For the 
minimizations, all parameters were fixed except Z300, the sur- 
face density at 300 AU. The adopted disk parameters are given 
in Table |2] Theses values come from the angularly reso lved CO 
interferome tric analysis performed bv lPietu et al.l (l2007h and for 
GO Tau by ISchaefer et al] ( |2009|) . The 3cr upper limits on the 
surface densities at 300 AU were then calculated from the for- 
mal error obtained from the best fit. The results are given in 
Table |3] The choice of the CO parameters to fit sulfur-bearing 
molecules is governed by the fact that CO results provide the 
best constraints on the gas disk structure and kinetic temperature. 
Moreover, the temperature derived from CO data can be consid- 
ered as representative of a significant fraction of the molecular 
layer so is well suited to determining upper limits on the surface 
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Table 3. Sulfur-bearing molecules: detections and 3cr upper lim- 
its. 
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Sources 




Sjoo (cm -) 








SO 


H2S 




cs 


DMTau 


< 7.5-10" 


< 1.4-10" 


3.5 


+ 0.1- 10'^ 


LkCalS 


< 1.9-10'2 


< 3.6-10" 


8.7 


± 1.6- 10'2 


MWC480 


< 2.5-10'^ 


<4.1-10" 


< 


8.4 -10" 


GOTau 


< 8.9 -10" 


< 1.8 -10" 


2.0 = 


t0.16- 10'- 



Notes. 



at 300 AU (modeled as I.{r) 
230o(''/300AU)"' '') derived from the 30-m data (except for 
3-2 in DM Tau) and the model DISKFIT. See text for details. 



Surface densities 



CS 



Table 4. Initial abundances of various species (/H) for the parent 
cloud. 



Element 


Abundance 


Ref 


Case* 


He 


0.09 


(1) 




N 


6.2(-5) 


(2) 




C+ 


1.7(-4) 


(2) 


A,C 




1.2(-4) 


(1) 


B,D 





1.4(-4) 


(3) 


A,C 




2.56(-4) 


(1) 


B,D 


s+ 


8(-8) 


(4) 


A,B 




8(-9) 




CD 


Si-^ 


8(-9) 


(4) 




Na+ 


2(-9) 


(4) 




Fe+ 


3(-9) 


(4) 




cr 


l(-9) 


(4) 






7(-9) 


(4) 




P + 


2(-10) 


(4) 





No tes. References : (PfWakelam & Herbsd (l20()8h. (2) | Jenkinsl j2009li . 
(3) lHincelin eFai] j2011h . and (4) lGraedel et all jl9823r * indicates to 
which chemical model the corresponding values apply. 



density of molecular species. Two parameters may affect £300: 
the exponent slope p and the outer radius Rout. 2^300 itself is rela- 
tively unsensitive to these parameters, but the derived error-bars 
can be significantly affected. On one hand, for DM Tau, using 
p - 0.5 instead of 1.5 decreases the error by a factor 4. On the 
other hand, further reducing Rout to 560 AU (the best fit value 
for CS in DM Tau) again increases the error by a factor 2. Thus 
our choice of values for Rout and p is rather conservative in view 
of the unknowns. 



3. Chemical analysis and modeling 

We used the NAUTILUS ( iHersant et al.ll2009h gas-grain model 
to study the sulfur chemistry in the four protoplanetary 
disks. NAUTILUS computes the abundance of 460 gas- 
phase and 195 surface species as a function of tim e us- 
ing the rate equation method (iHasegawa et alJ Il992h . The 
chemical network contains 4406 gas-phase reactions and 
1733 reactions involving grains (including adsorption and 
desorption processes and grain-surface reactions). The gas- 
phase network was updated according to the recommenda- 
tions from the KIDA experts (|http://kida.obs.u-bordeau xl .fr/l 
on November 2010). The full network can he found at 
|http://kida.obs.u-bordeauxl .fr/models/benchmark^OlO.dat We 
used silicate grains of 0.1 pm size. Th is small size agrees wit h 
the recent observational results from iGuilloteau et al.l (1201 ih . 



3 ^ 



2 r 



1 r 



10 15 20 25 
Temperature (K) 



30 




Fig. 2. Temperature, H2 density, and visual extinction vertical 
profiles in the four protoplanetary disks at 300 AU used for 
the chemical modeling. GO Tau and DM Tau have very simi- 
lar physical parameters so that we only show the curves for DM 
Tau. 



which show that in outer disks (r >70-100 AU) dust grains 
are dominated by small particles (ISM-like). This agrees with 
the theoretical results of !Birnst iel et al.l (1201 Ol) and suggests that 
large grains have akeady migrated to the inner disk regions 
(r <70 AU). Since small grains of O.lpm should remain dynam- 
ically coupled to the gas, this observational result also implies 
that the chemistry in the outer disk would be less affected by the 
dust settling than in the inner disk. For the ratio of extinction 
curve to the column density of hydrog en, we used the standard 
relatio n Ay = A^(//)/1.6-10-' following lWagenblast & Hartauisl 
Mo re details on the model can be found in lHersant et al.l 
(2009D and lSemenov et al.l (l2010l) . To obtain initial abundances 
for the chemistry of protoplanetary disks, we first compute the 
chemical composition of the parent molecular cloud. With this 
aim, we run NAUTILUS during 10^ yr for the following physi- 
cal conditions: gas and dust temperature are 10 K, H density is 
2x10'* cm"^ and Ay is 50 (models A and B). We checked higher 



densities (2 x 10 cm" ) and older clouds (10 yr), correspond- 
ing to models C and D and also tested smaller Ay (10). For the 
cloud, the species (listed in Table 4) are assumed to initially be 
in the atomic form except for hydrogen, which is entirely in H2, 
and we assumed four sets of elemental abundances. In models A 
and C, we adopted the oxygen elemental abundance required to 
reproduce the low O2 abundance observed in the cold ISM from 
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Table 5. Parameters used for the four models A, B, C, and D. 



Table 6. Model column densities 



Parameter 


Model A 


Model B 


Model C 


Model D 


Source 


Nso 


Nh,s 


Ncs 


C/O ratio 


1.2 


0.4 


1.2 


0.4 


MWC 480 


6.6-10'^ 


7.9-10'^ 


2.0-10'' 


Sulfur abundance 


8 X 10-** 


8 X 10-** 


8 X 10-' 


8 X 10-' 


LkCal5 


1.3-10'3 


6.8-10'3 


l.MO"' 


(/H) 










DM Tau " 


2.2-10'3 


6.2-10'3 


8.0-10'3 


Cloud H density 


2x 10-* 


2x 10** 


2x 10^ 


2x 10^ 


MWC 480 


7. MO" 


5.5-10'^ 


1.4-10'^ 


(H.cm-3) 










LkCal5 


1.3-10'^ 


9.M0'2 


l.MO'^ 


Cloud age 


105 


105 


10* 


10« 


DM Tau " 


2.4-10'2 


7.3-10'2 


8.9-10'2 



(yi's) 



iHincelin et al.l ( 1201 lb and the carbon abundance from iJenkinsI 
(I2009h for the ^Oph diffuse cloud, so that the C/O elemental ra- 
tio is 1.2. In models B an d D, we use the carbon and oxygen el- 
emental abundances from lWakelam & HerbstI (l2008l 1.2 x 10 "^ 
for C and 2.56 x lO ** for O), which gives a C/O ratio of 0.4. 
The choice of the sulfur elemental abundance "available" for the 
chemistry, i.e. not locked into the refractory part of the grains, 
is a widely studied problem (Ruffle et al. 199 9; Wakelam et aP 
I2OOI IScappini et al]l2003l: Ivan der Tak et"ani2003i) . No deple- 
tion of gas phase atomic sulfur is observed in diffuse clouds 
compared to the cosmic abundance of 10"^ (/H), in contrast to 
other elements. However, if such high elemental abundance is 
used to do the chemistry in dense clouds, one would produce 
too many S -bearing molecules, on which we have observat ional 
constraints dGraedel et al.lll982bl:IWakelam & Herbsi2008h . For 
this reason, depletion of sulfur in an unknown, more or less 
refractory form is assumed to happen very quickly when the 
molecular cloud is formed from the diffuse gas. To simulate this 
effect, the "free" elemental abundance of S is then decreased 
by an arbitrary amount. For simplicity, we have adopted the el- 
emental abundance from iGraedel et al.l (Il982bh in our Models 
A and B (which approximately reproduces the observations in 
dense clouds) and decreased it by a factor of 10 in models C and 
D. Table |5] also summarizes the parameters (cloud age, initial 
density, C/O ratio) for the four Models A, B, C, and D. We did 
numerous simulations to understand the effect of the parameters 
that we changed. In general, the S-bearing molecular abundances 
are not very sensitive to these parameters, and we present here 
the extreme cases. 

The abundances of chemical species in the gas and on the 
grains were used as initial composition for the chemistry of the 
disks, which is then integrated over 5x10^ yr, taken as the ap- 
proximate age of the studied disks. The computed column den- 
sities do not vary significantly after 10^ yr For the disk physical 
parameter s, we used the two layer parametric disk model de- 
sc ribed in Hersant et al.l (l2009t) . Its temperature profile is based 
on lDartois et al.l(l2003h . The dust temperature is equal to the gas 
temperature. The midplane is cold, with a temperature of 10 K 
for DM Tau, GO Tau, and LkCalS and 20 K for MWC480, con- 
sistent with the temperature derived from the '^CO 1-0 transi- 
tion. The molecular surface layer is warmer, as derived from 
'^CO 2-1 observations, with a temperature of 17.3 K for DM 
Tau and GO Tau, 12.7 K for LkCalS, and 27.7 K for MWC480 at 
300 AU. We assumed the same surface density for all four disks, 
with a value of 0.15 g cm^^ at 300 AU. The res ulting disk struc- 
ture is displayed Fig.l2l(see lHersant et al.ll2009l for details about 
the computation of the physical structure). For the UV radiative 
transfer in the disk, we consider only the vertical extinction (as- 
suming that only half of UV photons are scattered downwards 
by small grains located above our computing box). The UV flux 
given at 100 AU in Table[T]is decreasing as 1/r^ (where r is the 



Notes. Column densities (cm"-) of SO, CS, and H2S computed by the 
chemical model in the four protoplanetary disks at 300 AU from the 
central star for a disk 5 x 10'^ yr old. Top: Model A. Bottom: Model C. 
" Modeled abundances are the same for DM Tau and GO Tau. 



100 



z° 0.1 



0.01 



0.001 




Fig. 3. Ratio between the column densities of SO, CS, and H2S 
observed in the four protoplanetary disks at 300 AU from the 
central star and the values computed with NAUTILUS. For SO 
and H2S, they correspond to upper limits. The four Models (A, 
B, C, and D) correspond to the elemental abundances given in 
Tableland to the parameters given in Table|5] 



spherical radius) from the central star At a cylindrical radius of 
300 AU and a vertical boundary of four pressure scale heights 
(//), the UV scale factor becomes 33 for GO Tau, 17 for DM 
Tau, 120 for LkCalS, and 391 for MWC 480, expresse d in units 
of the Interstellar Radiation Field (ISRF. lDraineil978h . 

Using the chemical and physical parameters, we ran the four 
Models A,B,C and D in order to obtain the column densities of 
SO, H2S, and CS at 300 AU. The results for Models A and C 
are given in Table |6] and the comparison with the observations 
is shown in Fig. |3] where we present the ratio between predicted 
and observed column densities for the three molecules in the four 
sources for the four cases. 



4. Discussion 

In Fig. [3] except for CS, we only have upper limits on the ratios. 
DM Tau and GO Tau have similar model parameters so that the 
chemical models predict similar results for the two disks. We 
chose 5 X 10^ yr as the approximate age of the studied disks. 

Models A and B are not very sensitive to the density of 
the parent cloud and its age. The difference between H2S and 
CS column densities are small, except for MWC 480 with the 
largest UV flux. When the sulfur is less abundant by a factor of 
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Fig. 4. Modeled column densities for SO (diamond), H2S (cross) 
and CS (star) for the four sources in the case of Model C (C/O - 
1.2, initial cloud H density of 2 X 10^ cm"^ and age of 10^ yr). 
The observed column densities have been reported in grey (or 
color). Top panel: MWC480, medium panel: LkCalS, Bottom 
panel: GO Tau and DM Tau (same model). 



10 (Models C and D), the results are slightly more sensitive to 
the density and age of the cloud. 

In the case of Models C and D, corresponding to a cloud 
with an initial H density of 2 x 10^ cm"^ and an age of 10'' yr, 
the general agreement with the observations is better. The col- 
umn densities for CS agree with the observed values, and the 
ones obtained for SO agree with the upper limits. The column 



density for H2S, in contrast, is still too high by at least one order 
of magnitude. This result seems to indicate some missing pro- 
cesses for the destruction o f H2S on grains. H2 S is not formed 
efficiently in the gas-phase dHerbst et al.lll989l) . The gas phase 
H2S, in these simulations, comes from the non-thermal evapo- 
ration of solid H2S formed on grains by the hydrogenation of 
atomic sulfur (depleted from the gas phase). Its main destruc- 
tion path in the gas-phase is through the reaction with C^. At 
the high densities and relatively low temperatures encountered 
around disk mid planes (for H7S , the ev aporation temperature is 
around 40-50 K. lwakelam et alj (1201 lb ). H2S is Ukely to remain 
locked onto the grain surfaces and may react with other species 
preventing desorption of H2S. These grain surface reactions have 
not yet been inc orporated in chemica l models. In recent exper- 
imental studies. Rjarozzo et al. l (l2OT0h have shown that H2S on 
grains is easily destroyed by cosmic-ray particles, leading to the 
formation of C2S, SO2 and OCS on grains. In these experiments, 
most of the sulfur, however, may be in the form of a sulfur-rich 
residuum, which could be polyme rs of sulfur or amorph ous ag- 
gregates of sulfur, as suggested bv lWakelam et al.l ( |2004|) . 

We also checked the sensitivity to the C/O elemental ratio 
(C/0=1.2 for A and C and 0.4 for B and D). In Models B and 
D with C/O=0.4, we obtain an SO column density that is higher 
than the one of CS. The agreement with the observations is thus 
always less good than using the higher C/O value. 

We have superimposed in Fig|4]the observed column densi- 
ties on the predicted column densities for the three species in all 
sources for the Model C, which appears slightly better than the 
case of D. The column densities have been calculated at radii, 
50, 100, 200, and 300 AU. The three modeled molecular col- 
umn densities show the same radial trends in all sources with the 
exception of SO in MWC480 and DM Tau at 100 AU. The differ- 
ence in the SO abundances in the three sources can be explained 
by a balance between formation and destruction efficiencies of 
this molecule. At this radius, SO is formed by the neutral-neutral 
reaction S + OH — » SO + H and OH is formed by the associa- 
tion radiation O H- H — > OH + hv. Water molecules are much 
less abundant than H or O. It is the association between O and 
H that quantitatively produces OH. In MWC480 and LkCal5, 
the abundances of atomic O and H are larger than in DM Tau 
because of larger UV fluxes, producing more SO. In MWC480, 
however, SO is more efficiently destroyed than in LkCal5 be- 
cause of larger abundances of ionized atomic carbon lower in 
the disk atmosphere (C+ + SO ^ CO+ + S). 

We present in Fig|5]the vertical surface density profiles at ra- 
dius 300 AU for the sulfur-bearing molecules for the C model, 
but we observe the same trends with all models. We clearly see 
the location of the molecular layer at Z/// < 2, but the sulfur- 
bearing molecule abundances saturate at different heights. The 
SO vertical profile is sharp and peaks at a Z/H of 1.7 in DM 
Tau and GO Tau, 1.6 in LkCa 15 and 1.1 in MWC 480. The 
H2S column densities saturate slightly below in DM Tau, GO 
Tau, and LkCal5. In MWC 480, the profile is less peaked and 
spreads over 1.5 scale heights above the midplane. CS abun- 
dance has a rather broad maximum around two scale heights 
in all sources, except MWC 480, where it peaks at Z/H = 1.5 
with a narrow profile. The different behavior of MWC 480 can 
be understood because the UV flux is so high that all molecules 
dissociate above 1.7-2 scale heights, resulting in a narrower 
chemically rich layer and lower overall column densities. 

The relative abundanc es may be compared to those found 
in comets like Hale-Bopp (iBockelee-M orvan et al. 2004). In the 
coma of comets, the molecules are removed from the grains 
by desorption. Some species, like H2S, are directly desorbed 
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Table 7. Comparison at 300 AU between the observed ratios 
H2S/CO in the three disks of MWC480, DM Tau and LkCalS 
and the same ratios in comet Hale-Bopp 



Source 


^300 (cm"^) 


H2S/CO 


DM Tau 


3.7x10" 


< 4 X 10--' 


LkCalS 


1.8x10" 


< 10^^ 


MWC 480 


2.3x10" 


< 10--^ 


Hale-Bopp 




0.07 



Notes. The CO surface densities at 300 AU jPietu et al.ll2007h are ex- 
trapolated from '^CO and multiplied by 70. 



and most of them are short-lived because they are photodissoci- 
ated to give daughter molecules. Both CS and SO are produced 
by ph otodissociation of SO2 and CS2 dBockelee-Morvan et al.l 
I2OO4I) . Therefore, the quantity of H2S observed in Hale-Bopp 
is a direct tracer of the amount of H2S remaining currently on 
cometary grains, while the CS and SO observations most likely 
trace the reservoir on grains of CS2 and SO2. 

The ratios H2S/CO in Hale-Bopp and our disks are quoted 
in Table Q For our disk observations, we use the measured sur- 
face densities of C O derived from the interferometric analysis of 
iPietu et alJ (l2007h (in fact '^CO multiplied by 70, to avoid pos- 
sible opacity eff'ects). Our results are obtained at 300 AU (the 
region to which we are most sensitive) and trace the gas phase, 
while the Hale-Bopp measurements trace the current molecu- 
lar composition of the grain mantles probably formed around 5- 
30 AU in the protosolar nebula. The huge difference between 
observed gas ratio in the protoplanetary environment and the 
cometary ice ratio, at least a factor ~ 1000, suggests that grain 
surface chemistry on the comet may have actively modified the 
molecular content of the protoplanetary ice composition. 



5. Conclusions 

We have reported a new deep search for sulfur-bearing 
molecules (SO, H2S, and CS) in three TTauri disks (DM Tau, 
GO Tau, and LkCalS) and one Herbig Ae object (MWC480). 
CS was detected in the three TTauri disks but not in the Herbig 
Ae one. H2S and SO were not detected, but the upper limits are 
significant and allow us to make comparisons with models. For 
this purpose, we used the astro-chemistry code NAUTILUS. We 
find that better agreement with the observations is obtained for 
an initial cloud H density of 2 x 10^ c m"^ and age of 1 0^ yr with 
a C/O of 1.2 (as recently suggested bv lHincelin et alj ( 1201 lb ). 

Although it reproduces the SO and CS column densities rea- 
sonably well, our best model fails to reproduce the upper limits 
obtained on H2S by at least one order of magnitude, suggesting 
that a fraction of sulfur may be depleted in mantles or refractory 
grains. At the high densities and low temperatures encountered 
around disk midplanes, H2S may likely remain locked onto the 
grain surfaces, where it may react to form other species prevent- 
ing desorption of H2S. These kinds of grain surface reactions 
have not yet been incorporated in models. Our results empha- 
size the need for grain surface reactions in astrochemical models 
in presence of the high density and low temperature associated to 
a UV photon source, as is the case in protoplanetary disks. More 
sensitive observations with ALMA, are needed to provide astro- 
chemists with stronger constraints on the missing ingredients in 
disk chemical models. 
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